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ABSTRACT 

Context. The red supergiant (RSG) Betelgeuse is an irregular variable star. Convection may play an important role in understanding 
this variability. Interferometric observations can be interpreted using sophisticated simulations of stellar convection. 
Aims. We compare the visibility curves and closure phases obtained from our 3D simulation of RSG convection with C05B0LD to 
various interferometric observations of Betelgeuse from the optical to the H band in order to characterize and measure the convection 
pattern on this star. 

Methods. We use 3D radiative-hydrodynamics (RHD) simulation to compute intensity maps in different filters and we thus derive 
interferometric observables using the post-processing radiative transfer code 0PTIM3D. The synthetic visibility curves and closure 
phases are compared to observations. 

Results. We provide a robust detection of the granulation pattern on the surface of Betelgeuse in the optical and in the H band based 
on excellent fits to the observed visibility points and closure phases. Moreover, we determine that the Betelgeuse surface in the H 
band is covered by small to medium scale (5-15 mas) convection-related surface structures and a large (^30 mas) convective cell. In 
this spectral region, H2O molecules are the main absorbers and contribute to the small structures and to the position of the first null 
of the visibility curve (i.e. the apparent stellar radius). 

Key words, stars: betelgeuse - stars: atmospheres - hydrodynamics - radiative transfer - techniques: interferometric 

1, Introduction The position of Betelgeuse on the H-R diagram is highly 

uncertain, due to un certainty in its effective temperature. 
Betelgeuse is a red supergiant star (Betelgeuse, HD 39801, Ml- iLevesaue et"aD (l2005h used on e-dimensional MARCS models 



2Ia-Ibe) and is one of the brightest stars in the optical and near dOustafsson et al.1 120031 Il975h to fit the incredibly rich TiO 



1^ infrared. This star exhibits variations in integrated brightness, molecular bands in the optical region of the spectrum for sev- 

^ ^ surface features, and the depths, shapes, and Doppler shifts of eral RSGs. They found an effective temperature of 3650 K for 

; its spectral lines. There is a backlog of visual- wavelength obser- Betelgeuse. Despite the fact that they obtained a good agree- 

^ . vations of its brightness which covers almost a hundred years, ment with the evolutionary tracks, problems remain. There 

C3 I The irregular fluctuations of its light cur ve are clearly ape riodic is a mismatch in the IR colors that could be due to atmo- 

■ ■ ■ ' and rather resemble a series of outbursts. iKissetaD (|2006|) stud- spheric temperatu re inhomogeneities characteristic of convec- 

ied the variability of different red supergiant (RSG) stars includ- tion (iLevesque et a l. 2006). Also the distance of Betelgeuse has 

ing Betelgeuse and they found a strong noise component in the large uncertainities because of err ors related to the po sitional 

photometric variability, probably caused by the large convection movement of the stellar photocenter. lHarper et aP (l2008l) derived 

ceUs. In addition to this, the spectral line variations have been a distance of (197+45 pc) using high spatial resolution, mul- 

analyzed by several authors, who inferred the presence of large tiwavelength, VLA radio positions combined with Hipparcos 

granules and high convective velocities ( Josselin & Plez 2007; Catalogue Intermediate Astrometric Data. 
lGrayll200^ . iGrayl also found line bisectors that predominantly 

have reversed C-shapes, and line shape variations occurring at Betelgeuse is one of the best studied RSGs in term of multi- 

the 1 km s"Mevel that have no obvious connection to their shifts wavelength imaging because of its large luminosity and an- 

in wavelength. gular diameter. The existence of hot spots on its surface has 

been proposed to explain n umerous interferometric observa- 
tions with WHT and COAST (iBuscher et all 19901: 1 Wilson et alJ 
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|2004|) that have detected time- variable inhomogeneities in the 
brightness distribution. These authors fitted the visibiHty and 
closure phase data with a circular limb-darkened disk and 
zero to three spots. A large spot has been also detected by 
lUitenbroek et al.l (1 19981) with H ST. The non- spheri cal shape of 
Betelgeu se was also detected b y iTatebe et al.l (l2QQ7l) in the mid- 
infrared. iHaubois et al.l (l2009l) published a reconstructed image 
of Betelgeuse in the H band with two spots using the same 
data as presented in this work. iKervella et al.l (l2009l) resolved 
Betelgeuse using diff'raction-limited adaptive optics in the near- 
infrared and found an asymmetric envelope around the star 
with a bright plume extending in the southwestern region. They 
claimed the plume was either due to the presenc e of a convec- 
tive h ot spot or was caused by stellar rotation. lOhnaka et al.l 
(l2009h presented VLTI/AMBER observations of Betelgeuse at 
high spectral resolution and spatially resolved CO gas motions. 
They claimed that these motions were related to convective mo- 
tions in the upper atmosphere or to intermittent mass ejections 
in clumps or arcs. 

Radiation hydrodyn amics (RHP) simu lations of red super- 
giant stars are availab le (iFrevtag et al.l2QQ2|) to interpret past and 
future observations. (IChiavassa et al T 120091 hereafter Paper I) 
used these simulations to explore the impact of the granulation 
pattern on observed visibility curves and closure phases and de- 
tected a granulation pattern on Betelge use in the K band b y fit- 
ting the existing interferometric data of iPerrin et al.l (|2004|) . 

This paper is the second in the series aimed at exploring the 
convection in RSGs. The main purpose is to compare RHD sim- 
ulations to high-angular resolution observations of Betelgeuse 
covering a wide spectral range from the optical region to the 
near-infrared H band, in order to confirm the presence of con- 
vective cells on its surface. 



2. 3D radiation-hydrodynamics simulations and 
post-processing radiative transfer 

We employe d numerical simulations obtained using CO^BOLD 
(IFrevtag et al. 2002; Freytag 2003; Freyta g & Hofn er 2008) and 
in particular the model St35gm03n07 that has been deeply ana- 
lyzed in Paper I. The model has 12 M©, a numerical resolution of 
235^ grid points with a step of 8.6 7?©, an average luminosity over 
spherical shells and over time of L=93000±1300 L©, an eff'ective 
temperature of Teff =3490+13 K, a radius of 7?=832±0.7 7?©, and 
surface gravity log(g)=-0.337±0.001. This is our "best" RHD 
simulation so far because it has stellar paramet ers closest to 
Bet elgeuse {T^^ = 364 K. iLevesque et al.L 120051 and log(g)=- 
0.3. lHarpere"taD.l2008h . 

We used the 3D pure-LTE radiative transfer code 0PTIM3D 
described in Paper I to compute intensity maps from all the 
suitable snapshots of the 3D hydrodynamical simulation. The 
code takes into account the Doppler shifts caused by the convec- 
tive motions. The radiation transfer is calculated in detail using 
pre-tabulated extinction coefficients generated with the MARCS 
code dCrustafsson et al. 2008). These tables are functions of tem- 
perature, density and wavelength, and wer e computed with the 
solar composition of lAsplund et al.l (l2006h . The tables include 
the same extensive atomic and molecular data as the MARCS 
models. They were constructed with no micro-turbulence broad- 
ening and the temperature and density distributions are opti- 
mized to cover the values encountered in the outer layers of the 
RHD simulations. 



Table 1. Log of the observations. 



Date 


Telescope 


Filter (central A) 


October 7, 2005 


IOTA 


T/^XTT/^ 1 f\f\f\ 

IONIC - 16000 


ijcioDer o, zuuj 


TOT A 


ilJiML. - IDUUU 




IOTA 




October 11, 2005 


IOTA 


IONIC - 16000 


October 12, 2005 


IOTA 


IONIC - 16000 


October 16, 2005 


IOTA 


IONIC - 16000 


October 21, 1997 


COAST 


12900 A 


October 24, 1997 


COAST 


9050 A 


October 31, 1997 


COAST 


9050 A 


November 11, 1997 


COAST 


12900 A 


November 12, 1997 


COAST 


9050 A 


November 15, 1997 


WHT 


7000 A 


November 16, 1997 


WHT 


9050 A 


November 21, 1997 


COAST 


9050 A 


January 29, 2004 


COAST 


7500, 7820, 9050 A 


February 8, 2004 


COAST 


7500, 7820, 9050 A 


February 25, 2004 


COAST 


7500, 7820, 9050 A 


February 29, 2004 


COAST 


7500, 7820, 9050 A 


March 1, 2004 


COAST 


7500, 7820, 9050 A 


March 2, 2004 


COAST 


7500, 7820, 9050 A 



3. Observations 

The data presented in this work have been taken by two inde- 
pendent groups with different telescopes and they cover a large 
wavelength range from the optical to the near infrared. The log 
of the observations is reported in Table [T] 

3.1. Data at 16400 A. 

The H band data were acquired with the 3 tele scope inter- 
ferom eter IOTA (Infrared Optical Telescope Array, iTraub et aP 
I2OO3I) located at Mount Hopkins in Arizona. Light collected 
by three apertures (siderostats of 0.45m diameter) was spa- 
tially filtered by single mode fibers to clean the wavefronts, re- 
moving high frequency atmospheric corrugations that aff'ect the 
fringe contrast. Th e beams were then combined with IONIC 
(iBerger et al.ll2003l) . This integrated optics component combines 
3 input beams in a pairwise manner. Fringes were encoded in 
the time domain using piezo-electric path modulators, and de- 
tected with a near-in frared camera utilizing a PICNIC detector 
(iPedretti et al.ll2004l) . 

Betelgeuse was observed in the H band (16400 ± 1000 A, 
Fig.[B on 6 nights between 7th October and 16th October 2005. 
Five diff'erent configurations of the interferometer telescopes 
were used, in order to cover a large range of spatial frequencies 
between 12 and 95 arcsec"^ To calibrate the instrumental trans- 
fer function, observations of Betelgeuse were interleaved with 
observations of a reference (calibrator) star, HD 36167. 

Data reduction was carried out using an IDL pipeline 
(iMonnier et al ] l2004l:IZhaoetd]l2007h . In order to measure the 
closure phase, ^ ye took the phase of the complex triple prod- 
uct (bispectrum jBaldwin et al.Lri986h . The instrumental closure 
phase of I0NIC3 drifted less than 1 degree over many hours 
owing to the miniature dimensions of the integrated optics com- 
ponent. Both for the squared visibilities and the closure phase, 
the random errors were calculated with the bootstrap technique, 
in which a statistic is repeatedly re-estimated by Monte-Carlo 
sampling the original data with replacement. Full details of the 
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TiO 




Wavelength {f^m) 



GN CN,C0andH20 

2.5x10®^ r — ' — ' — ■ — ' — " — ' — ■ — ' — ' — ' — ■ — ' — " — ' — '~ 




Wavelength (/xm) 

Fig. 1. Transmission curves (red) of the bandpass filters used 
for Betelgeuse observations at COAST and WHT: central wave- 
lengths at 7000, 7500, 7820, 12900 A and nominal bandwidths 
(FWHM) of 100, 130, 50, 500, 1500 A respectively (lYoung et alJ 
l2004ll200Qh (the filter curve for the narrow bandpass centered at 
7820 A has been lost hence we assumed a top-hat filter); and for 
Betelg euse observations at IOTA: 16400 A with 1000 A band- 
width (iHaubois et al .1120091) . The black line is the synthetic spec- 
trum with its continuum computed from the RHD simulation 
snapshot of Fig.[8j the molecules that contribute the most in ev- 
ery filter have been highlighted. 



3.2.1. COAST data from 1997 

The COAST data were taken during October and November 
1997. Observations at 9050 A were made using the stan- 
dard beam-combiner and avalanche photodiode detectors 
[Baldwin et al ] (11994 . while 12900 A observations were ob- 
tained with_ase2arate_Q^ combiner optimized for JHK 
bands lYoung et al.l (Il998h . 

Observations of Betelgeuse were interleaved with observa- 
tions of calibrator stars, either unresolved or of small and known 
diameters. If at least three baselines were measurable and the at- 
mospheric coherence time was sufficiently long, closure phase 
measurements were also collected, by recording fringes on three 
baselines simultaneously. 

Data reduction was carried out using standard methods in 
which the power spectrum and bispec trum of the i nterfe rence 
fringes were averaged over each dataset [Burns et al.l (1 19971) . The 
resulting visibilities had formal fractional errors in the range 2- 
10 % of the values and the closure phases had typical uncertain- 
ties of 5-10°. Additional uncertainties of 10-20% were added 
to the visibility amplitudes to accommodate potential changes in 
the seeing conditions between observations of the science target 
and calibrator stars. 

3.2.2. WHT data from 1997 

The observations with the WHT used the non-redundant aper- 
ture masking method ([Baldwin et al.l 119861: iHaniff' etal.lll987h 
and employed a five-hole linear aperture mask. Filters centred at 
7000 A and 9050 A were used to select the observing waveband; 
only the 7000 A data are presented in this paper. The resulting in- 
terference fringes were imaged onto a CCD and one-dimensional 
fringe snapshots were recorded at 12-ms intervals. For each ori- 
entation of the r nask, the fringe data were reduced usin g stan- 
dard procedures (iHaniff' et al.lll987l: iBuscher et al.lll99Ql) to give 
estimates of the visibility amplitudes on all 10 interferometer 
baselines and of the closure phases on the 10 (linear) triangles 
of baselines. As for the COAST measurements, the uncertainties 
on the visibility amplitudes were dominated by calibration er- 
rors, which in this instance were unusually large (fractional error 
-30%). On the other hand, the calibrated closure phase measure- 
ments had typical errors of only 1-3°. The orientation and scale 
of the detector were determined by observations of two close 
visual binaries with well-determined orbits. 



observ ations and data reduction can be found in iHaubois et al.l 
(I2OO9I) . 

3.2. Data from 7000 to 12900 A. 

For this wavelength range, we used data tak en at two differen t 
epochs. The observations carried out in 1997 (lYoung et al .120001) 
were acquired with the Cambridge Optical Aperture Synthesis 
Telescope (COAST) on baselines up to 8.9 m (with central 
wavelengths/bandwidths of 9050/500 and 12900/1500 A) and 
by non-redundant aperture masking with the William Herschel 
Telescope (WHT) on baselines up to 3.7 m (with central wave- 
lengths/bandwidths of 7000/100 and 9050/500). The observa- 
tions carried out in 2004 were obtained with COAST on base- 
lines up to 6.1 m (bandpasses 7500/130, 7820/50 and 9050/500 
A). Fig.[T]shows all the filters used. 



3.2.3. COAST data from 2004 

The observations taken in 2004 (lYoung et al.l l2004l) were ac- 
quired with COAST using the standard beam combiner and fil- 
ters centered at 7500, 7820 and 9050 A with FWHM of 130, 50, 
and 500 Arespectively. The raw interference fringe data were re- 
duced using the same methods utilised for the 1997 COAST data 
to obtain a set of estimates of the visibility amplitude and closure 
phase for each observing waveband. 

4. Comparison of simulations and observations 

In this section, we compare the synthetic visibility curves and 
closure phases to the observations. For this purpose, we used 
all the snapshots from the RSG simulation to compute inten- 
sity maps with 0PTIM3D. These maps were normalized to the 

filter transmissions of Fig. [T] as: '^r^'^}^^^^ where 1(A) is the in- 
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tensity and r (A) is the optical transmission of the filter at a cer- 
tain wavelength. Then, for each intensity map, a discrete Fourier 
transform (FT) was calculated. The visibility V is defined as the 
modulus |z| of the complex Fourier transform z = x -\- iy (where 
X is real part of the complex number z and y its imaginary part) 
normalized to the modulus at the origin of the frequency plane 
Izol, with the phase defined as tan^ = The closure 

phase is defined as the phase of the triple product (or bispec- 
trum) of the complex visibilities on three baselines, which form 
a closed loop joining three stations A, B, and C. If the projection 
of the baseline AB is (wi,vi), that for BC is (w2,V2), and thus 
(ui + 1/2, vi + V2) for AC, the closure phase is: 

(f)c(uuVuU2,V2) = 
arg(y(t/i, vi) xV(u2, V2) x V*(ui + U2, vi + V2)). 

The projected baselines and stations are those of the observa- 
tions. 

Following the method explained in Paper I, we computed 
visibility curves and closure phases for 36 diff'erent rotation an- 
gles with a step of 5° from all the available intensity maps (^3.5 
years of stellar time), giving a total of ^^2000 synthetic visibili- 
ties and ^2000 synthetic closure phases per filter. 

4.1. Data at 16400 A 

We begin by comparing with the 16400 A data because this filter 
is centered where the H"^ continuous opacity minimum occurs. 
Consequently, the continuum-forming region is more visible and 
the granulation pattern is characterized by large scale granules 
of about 400-500 Rq (^60% of the stellar radius) evolving on 
a timescale of years (Fig. 4 in Paper I). On the top of these 
cells, there are short-lived (a few months to one year) small- 
scale (about 50-100 R©) structures. The resulting granulation 
pattern causes significant fluctuations of the visibility curves and 
the signal to be expected in the second, third and fourth lobes 
deviates greatly from that predicted by uniform disk (UD) and 
limb-darkened disk (LD) models (Fig. 11 in Paper I). Also the 
closure phases show large departures from and ±;r, the values 
which would indicate a point- symmetric brightness distribution. 

Within the large number of computed visibilities and closure 
phases for this filter, we found that some match the observation 
data very well (Fig.[2l). We selected the best-fitting snapshot min- 
imizing the function: 



^ N 



-M 



0C/ 



(1) 



where Vi is the observed visibility amplitude data with its cor- 
responding error cry.. My. is the synthetic visibility amplitude 
at the same spatial frequency, (pd is the observed closure phase 
with corresponding error cr^^ , and M^^ the synthetic closure 
phase for the observed UV coordinates. The best matching vis- 
ibilities and closure phases correspond to a particular snapshot 
and rotation angle. 

In Fig.[3]the simulation has been scaled to an apparent diam- 
eter of -45.1 mas in order to fit the data points in the first lobe, 
corresponding to a distance of 172.1 pc for the simulated star. 
The angular diameter is slightly larger t han the limb-darkene d 
diameter of 44.28+0.15 mas found bv iHaubois et al.l (l2009l) . 
Our distance is also in agreement with Harp er et al .l (l2008hT who 
reported a distance of 197 ± 45 pc. Using Harper et al.'s distanc e 
and an eff'ective temperature of 3650 K (iLevesque et al.l [2005b . 
the radius is R = 890 ± 2007?©, neglecting any uncertainty in 
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Fig. 2. Top and central panels: the best-matching synthetic 
squared vis ibility (grey) comp ared to the observations of 
Betelgeuse (iHaubois et al.ll2009l red). The reduced is indi- 
cated. The dotted lines are the maximum and minimum am- 
plitude of the visibility fluctuations as the snapshot is rotated. 
Bottom panel: closure phase (in degrees) versus the maximum 
projected baseline length for each baseline triplet. The observa- 
tions (red dots with error bars) are fitted with the same intensity 
map and rotation angle used for the visibilities (black crosses). 
The axisymmetric case is represented by the solid black lines. 
The dotted lines are the maximum and minimum of the closure 
phase variation as the snapshot is rotated. 
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Teff . On the other hand, using H arper et al. 2008 d istance and 
the apparent diameter of 45 mas (iPerrin et al ] |2004 . the radius 
is R = 950 + 200/?Q. All thes e results match evolutionary tracks 
by iMeynet & Maede3 (l2003l) for an initial mass of between 15 
and 25 M©. The radius (R ^ 8327?©, see Section ^ and the 
effective temperature (Tq^ ^ 3490 K) of our 3D simulation 
are smaller because the simulations start with an initial model 
that has a guessed radius, a certain envelope mass, a certain 
potential profile, and a prescribed luminosity. However, during 
the run the internal structure relaxes so something not to far 
away from the initial guess (otherwise the numerical grid is 
inappropriate). The average final radius is determined once 
the simulation is finished. Therefore, since the radius (and the 
effective temperature) cannot be tuned, the model is placed at 
some distance in order to provide the angular diameter that 
best matches the observations. Finally, within the error bars 
our model radius agre es with all other da ta derived using the 
distance determined by lHarper et al.l (l2008l) . 

Our RHD simulation provides a better fit than u niform disk 
and limb-darkened models used by iHaubois et al.l in all lobes 
of the visibility function. The departure from circular symmetry 
is more evident at high spatial frequencies (e.g., the fourth lobe) 
where the visibility predicted from the parametric model is lower 
than the observed data. The small-scale convection-related sur- 
face structures are the cause of this departure and can only be 
explained by RHD simulations that are permeated with irregular 
convection-related structures of different size. 



Filter 16400 ±1000 A 




-20 -10 10 SO 

X [mas] 



Fig. 3. Map of the linear intensity in the IONIC filter. The range 
is [0;3.1 X 10^] ergcm"^ s"^ A"^ The stellar parameters of this 
snapshot are: L = 92 300 L©, R = 832.47?©, T^^ = 3486 K and 
log(g) = -0.34. The simulation has been scaled to an apparent 
diameter of -45.1 mas, at a distance of 172.1 pc in order to fit 
the data points in the first lobe. 

Also the closure phases display a good agreement with the 
simulation indicating that a possible solution to the distribution 
of the inhomogeneities on the surface of Betelgeuse is the in- 
tensity map of Fig. [3] (though the reconstructed images found by 
IHaubois et afl (120091) are more probable). 



This is the first robust confirmation of the physical origin of 
surface granulation for Betelgeuse, follo wing on from t he de- 
tection in the K band (Paper I). Recentlv. IHaubois et all (l2009l) 
were able to reconstruct two images of Betelgeuse, using the 
data presented in this work, with two different image recon- 
struction algorithms. The i mage reconstructed with WISARD 
(Meimon et aDl2009L 120051) is display ed in Fig. H (left). Both 
reconstructed images in IHaubois et al.l paper have two spots of 
unequal brightness located at roughly the same positions near 
the center of the stellar disk. One of these spots is half the stellar 
radius in size. Fig. |4] shows a comparison of the reconstructed 
image to our best fitting snapshot of Fig. O Fainter structures 
are visible in the synthetic image (right panel) while the recon- 
structed image (left panel) is dominated by two bright spots. 
Moreover, the bigger spot visible in the reconstructed image is 
not present in our synthetic image, whereas there is good agree- 
ment in term of location with the smaller spot located close to 
the center. However, it is possible that the synthetic map does 
not match exactly the location of the spots because it cannot per- 
fectly reproduce the closure phase data. 

4.1 .1 . Molecular contribution to the visibility curves 

It is important to determine what molecular species contribute 
the most to the intensity absorption in the stellar atmosphere. For 
this purpose, we used the best fitting snapshot of Fig. [3] and re- 
computed the intensity maps in the IONIC filter using only CO, 
CN and H2O molecules, because they are the largest absorbers 
at these wavelengths (Fig. 1 of this work and Fig. 3 of Paper I). 
The intensity maps displayed in Fig. \5\ of this paper (top row) 
should be compared to the original one in Fig. [3] which accounts 
for all the molecular and atomic lines. The surfaces of CO and 
CN maps clearly show the granulation pattern and they are spot- 
free. However, the H2O map shows dark spots which can be also 
identified in the original intensity map. 

We also calculated visibility curves from these molecular in- 
tensity maps using the same rotation angle used to generate the 
synthetic data in Fig. [2l Figure \5\ (bottom row) shows that the 
H2O visibility is the closest to the original one both at low and 
high spatial frequencies. We conclude that: (i) in the first lobe, 
the H2O visibility is smaller than the CO and CN visibilities. 
Thus the radius of the star is dependent on the H2O contribu- 
tion, (ii) At higher frequencies, only the H2O visibility can fit 
the observed data whereas the CO and CN visibilities fit poorly. 

4.1 .2. Size distribution on the stellar surface 

We also characterized the typical size distribution on the stel- 
lar surface using interferometric observables. The large range 
of spatial frequencies of the observation (between 12 to 95 
arcsec"^) is very well suited for this purpose. 

Out aim is to visualize the energy within the signal as a func- 
tion of spatial frequency. After the computation of the Fourier 
Transform, FT, we obtain I(u,v) = FT [I(x,y)] where I(u,v) 
is the best matching intensity map of Fig.O The resulting com- 
plex number / (w, v) is multiplied by low-pass and high-pass fil- 
ters to extract the information from different spatial frequency 
ranges (corresponding to the visibility lobes). Finally an inverse 
Fourier Transform, FT, is used to obtain the filtered image: 
/filtered (x,y) = FT [/(w, v) • filter]. 

Fig. [6] (top left panel) shows the filtered image at spatial fre- 
quencies, y, corresponding to the first lobe. Due to the fact that 
we cut off the signal at high spatial frequencies, the image ap- 
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20 -10 10 
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Fig. 4. Left panel: recontructed image from iHaubois et alJ (l2QQ9l) . Right panel: our best fitting 3D simulation snapshot of Fig. [3] 
convolved with a 6x6 mas PSF derived by fitting the central peak of the interferometric dirty beam. The intensities in both panels 
are normalized to the range [0, 1] and some contour lines are indicated (0.2, 0.5, 0.7, 0.8, 0.9 of the peak brightness). 
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pears blurry and seems to contain only the information about the 
stellar radius. However, the top right panel displays the signal 
related to all the frequencies larger than the first lobe: in this 
image we clearly miss the central convective cell of ^ 30 mas 
size (60% of the stellar radius) visible in Fig. O Thus, the first 



lobe also carries information on the presence of large convective 
cells. 

Fig. [6] (bottom row) shows the second lobe with convection- 
related structures of ^ 10-15 mas, (30% of the stellar radius), and 
the third and fourth lobes with structures smaller than 10 mas. 
We conclude that we can estimate the presence of convection- 
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related structures of different size using visibility measurements 
at the appropriate spatial frequencies. However, only imaging 
can definitively characterize the size of granules. A first step 
in this direction has been carried out in Berger et al. 2010 (to 
be sumbitted soon), where the image reconstruction algorythms 
have been tested using intensity maps from this RHD simula- 
tion. In the case of Betelgeuse, we have fitted its interferometric 
observables between 12 and 95 arcsec"^ and thus inferred the 
presence of small to medium scale granules (5 to 15 mas) and a 
large convective cell (^^ 30 mas). 

4.2. Data from 7000 to 12900 A 

The simulated surface in the optical to near-infrared region 
shows a spectacular pattern characterized by dark spots and 
bright areas. The brightest areas can be up to 50 times more 
intense than the dark ones. In addition, this pattern changes 
strongly with time and has a lifetime of a few weeks. In the 
wavelength region below ^ 1 yum, the resulting surface pat- 
tern, though related to the granulation below, is also connected 
to dynamical eff'ects. In fact, the light comes from higher up in 
the atmosphere where the optical depth is smaller than one and 
where waves and shocks start to dominate. In addition to this, 
the emerging intensity depends on the opacity run through the 
atmosphere and veiling by TiO molecules is very strong at these 
wavelengths (Fig.[T]). 

The filters centered in the optical part of this wavelength 
range are characterized by the presence of strong molecular lines 
while the infrared filter probes layers closer to the continuum- 
forming region. Observations at wavelengths in a molecular 
band and in the continuum probe diff'erent atmospheric depths, 
and thus layers at diff'erent temperatures. They provide impor- 
tant information on the wavelength-dependence of the limb- 
darkening and strong tests of our simulations. 

Since the observations have been made at two diff'erent 
epochs, we fitted th em individually: (i) the data t aken in 1997 
('Young et all l2000h . and (ii) the data from 2004 (lYoung et al.l 
2004). We proceeded as in Section |4Jl Note that the filter curve 
for the 7820 A bandpass was lost and so a top-hat was assumed 
instead; we tested the validity of this assumption by replacing 
the known optical-region filter curves with top-hats, which did 
not aff'ect the synthetic visibility and closure phase data signifi- 
cantly. Within the large number of computed visibilities and clo- 
sure phases for each filter, we found that there are two snapshots 
of the simulated star, one for each epoch, that fit the observa- 
tions. At each epoch the same rotation angle of the snapshot was 
found to fit all of the observed wavebands. 

Fig.[7]displays the comparison to the data taken in 1997. The 
7000 A synthetic image corresponds to a region with strong TiO 
absorption (transition A^O-X^A(7), see Fig.[T]). This is also true 
for the 9050 A image (transition E^H - X^A(e)) but in this case 
the TiO band is less strong. The relative intensity of TiO bands 
depend on the temperature gradient of the model and change 
smoothly from one snapshot to another. The map at 12900 A is 
TiO free and there are mostly CN lines: in this case, the surface 
intensity contrast is less strong than in the TiO bands. 

Fig. [8] shows the comparison to the data taken in 2004. 
There, the filters used span wavelength regions corresponding to 
TiO absorption bands with diff'erent strengths centered at 7500A 
7800A and 9050A (transitions A^^-X^A(y) and E^n-X^A(6), 
see Fig. [T]). Again, the same snapshot fitted the whole dataset 
from the same epoch. 



The departure from circular symmetry is more evident than 
in the H band, with the first and the second lobes already show- 
ing large visibility fluctuations. The RHD simulation shows an 
excellent agreement with the data both in the visibility curves 
and closure phases. 

However, within the same observation epoch we had to scale 
the size of the simulated star to a difl'erent apparent diameter 
at each observed wavelength. For example, in Fig. [8] the appar- 
ent diameter varies from 47.3 to 52.6 mas. Thus we infer that 
our RHD simulation fails to reproduce the TiO molecular band 
strengths probed by the three filters. As already pointed out in 
Paper I, our RHD simulations are constrained by execution time 
and they use a grey approximation for the radiative transfer that 
is well justified in the stellar interior and is a crude approxima- 
tion in the optically thin layers; as a consequence the thermal 
gradient is too s hallow and weakens th e contrast between strong 
and weak lines (IChiavassa et al.|[2006b . The intensity maps look 
too sharp with respect to the observations. The implementa- 
tion of non-grey opacities with five wavelength bins employed 
to describe the wavelength dependence of radiation fields (see 
Ludwig et al. 1994; Nordlund 1982, for details) should change 
the mean temperature structure and the temperature fluctuations. 
The mean thermal gradient in the outer layers, where TiO ab- 
sorption has a large efl'ect, should increase. Moreover, in a next 
step, the inclusion of the radiation pressure in the simulations 
should lead to a difl'erent density/pressure structure with a less 
steep decline of density with radius. We expect that the intensity 
maps probing TiO bands with difl'erent strengths will eventually 
show larger diameter variations due to the molecular absorption 
as a re sult of these r efinem ents. 

In I Young etaP (l2000l) . the authors managed to model the 
data in the 7000 A filter with two best-fitting parametric mod- 
els, consisting of a circular disk with superimposed bright fea- 
tures (Fig. [21 central panel) or dark features (left panel). We 
compared these parametric models to our best-fitting synthetic 
image of Fig. [7] (top row). The convolved image, displayed in 
Fig. [9] (right panel), shows a better qualitative agreement with the 
bright features parametric model. In fact, 3D simulations show 
that the surface contrast is enhanced by the presence of signif- 
icant molecular absorbers like TiO which contribute in layers 
where waves and shocks start to dominate. The location of bright 
spots is then a consequence of the underlying activity. 



5. Conclusions 

We have used radiation hydrodynamics simulations of red super- 
giant stars to explain interferometric observations of Betelgeuse 
from the optical to the infrared region. 

The picture of the surface of Betelgeuse resulting from this 
work is the following: (i) a granulation pattern is undoubt- 
edly present on the surface and the convection-related struc- 
tures have a strong signature in the visibility curve and clo- 
sure phases at high spatial frequencies in the H band and on 
the first and second lobes in the optical region, (ii) In the 
H band, Betelgeuse is characterized by a granulation pattern 
that is composed of convection-related structures of diff'erent 
sizes. There are small to medium scale granules (5-15 mas) 
and a large convective cell (^^30 mas). This supports previ- 
ous detections carried out with our RHD simulation in the K 
ban d (Paper I) using parametric models and the same dataset 
( Haubois et al l 120091). kervella et al. ( 2009) with VLT/NACO 
observations and lOhnaka et al.r(l2009l) with VLTI/AMBER ob- 
servations. Moreover, we have demonstrated that H2O molecules 
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Fig. 6. Intensity maps filtered at diff'erent spatial frequencies corresponding to the lobes of the visibility curve shown in Fig. [21 The 
images are normalized between [0,1]. 



contribute more than CO and CN to the position of the vis- 
ibility curve's first null (and thus to the measured stellar ra- 
dius) and to the small scale surface structures, (iii) In the op- 
tical, Betelgeuse's surface appears more complex with areas up 
to 50 times brighter than the dark ones. This picture is the conse- 
quence of the underlaying activity characterized by interactions 
between shock waves and non-radial pulsations in layers where 
there are strong TiO molecular bands. 

These observations provide a wealth of information about 
both the stars and our RHD models. The comparison with the ob- 
servations in the TiO bands allowed us to suggest which approx- 
imations must be replaced with more realistic treatments in the 
simulations. New models with wavelength resolution (i.e., non- 
grey opacities) are in progress and they will be tested against 
these observations. From the observational point of view, further 
multi-epoch observations, both in the optical and in the infrared, 
are needed to assess the time variability of convection. 
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Fig. 9. Left and central panels: Brightness distributions of the best fitting parametric models from I Young etaP (|200Q|) with dark 
features (left) and bright features (center). Right panel: our best fitting 3D simulation snapshot of Fig. [7] (row marked with (A)) 
convolved with a 10x10 mas PSF (i.e. the size of the bright features in the parametric model). The intensities in both panels are 
normalized to the range [0, 1] and some contour lines are indicated (0.1,0.2,0.5,0.7). 



